First principles study of iron-based molecule grafted on graphene 
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Motivated by recent experimental studies on single molecular magnets grafted on graphene and 
single walled carbon nanotubes, we investigate the structural, electronic, and magnetic properties 
of an iron based magnetic molecule grafted on a graphene sheet using ah initio calculations. In 
particular, the induced charge transfer and magnetization are described in terms of the coupling 
between the molecule and the graphene orbitals. This interaction and its effects on graphene elec- 
tronic properties are determined and discussed in view of the potential utilization of graphene in 
spintronics. 



I. INTRODUCTION 



Graphene-based systems have recently attracted much 
interest from both experimental and theoretical as- 
pects. Its chemical inertness, hydrophobic behavior, 
large electron mobility, scalable production and intrin- 
sically low spin-orbit coupling makes it a very promi- 
sing candidate for sensor^, electronicJ^, spintronicPSi 
and nanomechanic^^. Experimentally, the graphene ca- 
pability to detect small charge transfer effects by graf- 
ting molecule^i^, but also to induce superconductivitj!^, 
as well as Kondo effeclP^ by deposing metals have been 
investigated. More recently, surface-enhanced Raman si- 
gnal has been r ealize d on pyrene-based molecules graf- 
ted on graphen ci^^ l ^^ l, allowing a detection down to the 
limit of few isolated molecules. Detection of single mo- 
lecule magnelP'^ by using graphene or carbon nano- 
tubes could also be a way to probe magnetic properties at 
the single molecule level. Theoretically, ah initio calcula- 
tions on magnetic pr opert ies of graphene based materials 
have been performecpm^l. In particular, ab initio studies 
of adatoms on graphen^i^HlSI magnetic Co dimer^SI and 
even small magnetic molecules such as O2, NO or NOjp' 
have been reported. 

The interaction between the graphene monolayer and a 
magnetic molecule is an important point in order to have 
insights on local effects such as induced magnetization 
and charge transfer, which are of interest for spintronic 
and nano-electronic applications. In the present work, we 
chose the iron tetraphtalic acid molecule - FeTPA^^ - 
deposited on monolayer graphene as a case study of this 
interaction. Indeed, Gambardella et alW^ have succeeded 
in the manipulation of the magnetic anisotropy of a su- 
pramolecular assembly of FeTPA^ self assembled on a 
Cu surface. However, to speed up our DFT calculations, 
we modeled FeTPA^ by explicitly studying its magnetic 
core Fe{OH)i, as we will discuss in the following, where 
this choice will be justified on the basis of structural, 
electronic, and magnetic properties. 



II. COMPUTATIONAL DETAILS 

First principles DFT calculations were performed 
within the plane wave approximation (PW) and pseu- 
dopotential sche me, a s implemented in the quantum 
ESPRESSO cod£2i!22! We adopt a Perdew-Burke- 
Ernzerhof (PBE)^'' gradient corrected functional for 
the exchange and correlation potential, and ultrasoft 
pseudopotential technique is used to describe C, H, O 
and Fe atoms^"^. Electronic wave functions are written 
in terms of plane waves, with an energy up to 30 Ry, 
which is sufficient to ensure convergence of structural, 
electronic and magnetic properties. The electronic occu- 
pation is computed using Fermi-Dirac distribution, with 
a smearing parameter of 136 meV, which corresponds 
to an electronic temperature of 1578 K. The Brillouin 
zone integration is performed with a uniform k points 
grid of (15x15x1). The convergence of the electronic 
and spin properties with the electronic temperature 
and k-points grid has been tested. Our choice of 1578 
K with (15x15x1) k-points grid is within the range 
of typical electronic temperatures used to describe 
graphitic systems, and it is motivated by the fact that 
such values ensure converged results with reasonable 
computational effort. Structures are relaxed until forces 
are below 0.02 eV/A. We used a supercell composed by 
(6x6) repetition of a graphene unitary cell in order to 
ensure a negligible interaction between periodic images 
of the adsorbed molecule. Spin polarized calculations are 
also performed, and we consider polarization along the z 
axis (i.e. perpendicular to graphene plane). 



A. Modeling FeTPAi through a smaller molecule 

The FeTPA^ molecule is difficult to simulate due to 
its large number of atoms (33 atoms) and, especially, of 
the size of the graphene supercell that should contain 
it. As mentioned above, however, the Fe{OH)4 molecule 
mimics well the Fe{TPA)4 molecule. This smaller, case- 
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Figure 1. (Top and middle panels) Structural description of : 
a) Fe{TPA)i isolated molecule, b) Fe{OH)4, isolated mole- 
cule, c) hybrid system composed of graphene and Fe{OH)4. 
Iron and oxygen atoms are indicated in red, hydrogen atoms 
in blue, and carbon atoms in yellow. The high symmetry ad- 
sorption sites are indicated with respect to the hexagonal cell 
of graphene, (d). Bottom panel : charge density map within 
xy plan for e) Fe{TPA)4 and f) Fe{OH)i. 



study system contains the same FeO^ core as the bigger 
one, but the four large pyrenic arms are simply replaced 
by hydrogen atoms. As shown in the top panels of Fig. 
[T] the two molecules have a similar planar structure. The 
charge distribution around the iron core are computed 
projecting all the electronic wave functions tpi^a- on the 
atomic orbitals centered in each given atom through the 
equation : 
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Table I. Calculated Lowdin charges (see eq[T]) and their de- 
composition onto atomic-like orbitals, for Fe and O atoms 
in Fe{OH)i and Fe(TPA)4 molecules. O* stands for first 
neighbors of Fe atom. Polarization is computed along the z 
direction for Fe and O atoms as the difference between spin 
up and down charges. 



down charges. As shown in Table |Tj the spin polarization 
of Fe atoms is very similar, being 3.20 and 3.10 /is in 
Fe{OH)i and Fe{TPA)i molecules respectively. In both 
cases this polarization comes from unpaired electron in 
3d iron orbitals. The polarization of the O atoms, which 
comes from 2p orbitals, is slightly higher in the Fe{OH)4 
molecule. The total polarization of the core is 3.94 and 
3.59 fiB in Fe[OH)4^ and Fe(TPA)4, respectively, which 
corresponds to a variation of around 10%. We can thus 
safely consider these two molecules equivalent from the 
spin magnetic point of view in the void ; given the relative 
inertness of graphene sheets with respect to full metallic 
surfaces, we can reasonably extrapolate a very similar 
behavior when grafted on a graphene sheet. 



E 



(1) 



where (f)^^ is the atomic- like orbital for atom /3, labelled 
with its energy and angular momentum quantum number 

TO, I. 

We observe in Table |l] that the charge associated to 
the core iron atom (Lowdin charge) remains very similar 
in both cases. The only slight charge difference between 
the two cases concerns the O atoms in Fe{OH)i having 
around 0.2 electrons more than in Fe{TPA)i. We also 
computed the spin polarization along the z direction for 
Fe and O atoms, as the difference between the spin up and 



B. Calculation details of the hybrid system 



We use, to study the hybrid system composed by the 
Fe{OH)i molecule adsorbed on graphene, an hexagonal 
supercell composed of 36 unit cells (FigjlJ;) with a lattice 
parameter in the xy plane of ao = 14.76 A, and of 11.10 A 
in z direction. Due to Brillouin zone refolding, the Dirac 
(K) point of graphene refolds on the F point. The three 
high symmetry adsorption sites of FeiOH)^ on graphene 
are shown in Fig{lji. We find that the most energetically 
favorable adsorption site is the top site, and only this one 
will then be considered in the following. 
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III. RESULTS AND DISCUSSION 



A. Structural properties 
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Figure 2. (a) Binding energy of the hybrid system as a func- 
tion of the distance along z between the graphene surface and 
the iron atom of the Fe(0H)4 molecule, (b) Displacements 
(Az*^*) of graphene C atoms positions along z with respect 
to the isolated case, as a function of distance, within the xy 
plane, to the adsorption (top site) Ctop atom. The horizontal 
arrow represents the planar extension of Fe{OH)i molecule. 
Two fixed values of zfc-g are considered : the equilibrium 
distance 4.1 A (filled circles) and the characteristic distance 
3.1 A (empty triangles) 



We performed the structural relaxation of the hy- 
brid system by setting a fixed distance zpe-G between 
the Fe atom in the FeOHi molecule and the graphene 
plane. According to our calculations, top adsorption site 
is found to be the most stable from an energetic point 
of view and we defined ipe-G as the distance along z 



between the Fe atom and the first C atom below, which 
we indicate as Ctop- We calculated the binding energy 
as a function of this distance (see Fig. [2]), as Ehinding — 

Egraphene+FeOHi — Egraphene — EpeOHi- The obtained 

equilibrium distance for the system is z^g_g,=4.1 A, and 
the corresponding binding energy is Ehinding — —0.36 
eV. However, it is well known that ab initio DFT cal- 
culations do not take into account Van der Waals inter- 
actions, which are of course important in weakly inter- 
acting systems. As a result, the distance we obtain is li- 
kely to be overestimated. A properly ab initio correction 
should explicitly take into account these contributions in 
the DFT functional. However, this procedure is rather 
cumbersome and computationally costly, and a common 
approximation consists in carrying out the calculations 
at the experimental equilibrium distance. The most fla- 
grant case is graphite, a 3D crystal formed by layers of 
graphene, whose experimental interplane equilibrium dis- 
tance is about 3.35 A, while the DFT equilibrium value 
is more than 5 A. In an analogous way, in the following 
we always compare results obtained for the DFT equi- 
librium distance Zp*g_g,=4.1 A with results obtained at 

a shorter distance, zpe-c— 3.1 A which is comparable 
to the characteristic interaction distance observed in sp^ 
carbon materials^^ 

In both cases, the flatness of the graphene plane is slightly 
perturbed by the presence of the molecule. While C-C 
bond distances are practically unaffected (< 0.1 %), the 
z coordinate of C atoms in the relaxed hybrid system 
changes by as much as 1 %, as shown in Fig. [l] This sur- 
face modulation extends up to a distance of around 5.5 
A whether zp^^Q equals 4.1 or 3.1 A. 



B. Electronic properties 

The charge transfer between graphene layer and the 
molecule can be quantified by computing the difference 
between the total charge of the graphene sheet in the hy- 
brid system and of isolated graphene. We obtain a charge 
difference of -0.63 electrons, revealing an induced hole 
doping on graphene. This charge transfer is equal to a 
surface charge density of about -3.4 10^'^ cm~^, which is 
coherent with common predictions and experiments in- 
volving graphend^MSSl 

We also evaluate the local electronic charge transfer for 
each C atoms in the graphene plane and extract contri- 
butions from the different orbitals by projecting the elec- 
tronic eigenstates on the atomic wavefunctions as in Eq. 
[l] At a molecule-graphene distance of 4.1 A we observe 
that the charge transfer from the graphene sheet origi- 
nates mostly (98%) from the 2p orbitals of C atoms. A 
maximum charge transfer is observed for 2p states of C 
atoms at 2.5 A from the Fe center, while at distances grea- 
ter than 6 A the spatial charge transfer remains constant. 
For zpe-G — 3.1A, a similar behavior is observed. From 
this result we determine an interaction length of about 6 
A, which is in agreement with the one estimated by the 
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analysis of the surface corrugation (see Fig[T]). Moreover, 
the comparison between the electronic band structure of 
isolated graphene and of the hybrid system (not repor- 
ted) indicates that, despite a shift of the Fermi level, the 
perturbation induced by the molecule is negligible, as also 
shown by the non-dispersive character of the molecular 
states. 

This analysis of the charge transfer and of the electronic 
band structure of the hybrid system suggests that the 
quasi-metallic character of graphene is not perturbed by 
its interaction with the grafted magnetic molecule. This 
result is promising in view of the use of graphene as a 
sensitive detector. 



C. Magnetic properties 

As mentioned above, we also performed spin polarized 
calculations on isolated molecule and on the hybrid 
system. The obtained magnetic moment for the isolated 
molecule Fe{0H)4 is about 3.94 /is, whereas in the 
hybrid system it is about 3.34 fi^. The magnetization 
difference between the two systems is about 0.6 /xg indi- 
cating that the interaction with the graphene substrate 
generates 0.6 "paired" electrons in this hybrid system 
which is coherent with the estimated charge transfer of 
0.63 electrons found in the previous section. We then 
analyze this variation of the molecular magnetization 
within the hybrid system spin states occupations. Thus, 
we calculated the projected density of states (PDOS) for 
spin up and spin down states within the hybrid system. 
More precisely, we focused on the projection into 2pz 
states of carbon atoms, which are the most sensitive to 
their environment. 
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Figure 3. (a) Difference between PDOS (see text) for spin 
up and down states witfiin tfie liybrid system. Projections are 
performed onto carbon 2pz atomic orbitals. Carbon atoms 
are labelled by their neighbouring rank, in the xy plane, with 
respect to Ctop. The sketch represents the C atoms of the 
graphene. Fe{OH)4, is located at 3.1 A vertically above the 
black Ctop atom and the different colors indicates the different 
neighboring rank. PDOS2p^ (t) — PDOS2p, {i) is averaged for 
every C* atom at the same distance from Ctop- (b) Spatial 
distribution of the magnetization An(r) = n{f)^ ~ n(r)j_ wi- 
thin hybrid system for two distances zj?e-G = 4.1 A (top) and 
zpe-G ~ 3.1 A (bottom). 



In Fig. [3^ we report the difference between the PDOS 
for spin up and down states for different C atoms, la- 
belled with their distance in the xy plane with respect 
to Fe atom in the Fe{OH)4 molecule. This figure shows 
an energy range [-3 ; -1.5] eV where the PDOS for spin 
up and spin down differs and, in particular, for carbon 
atoms located up to the S*'' neighbour (ie. 5 A from Ctop 
atom). This spin- interaction distance is consistent with 
the previously discussed structural and electronic inter- 
action lengths of the hybrid system. 
To go further, we analyze the mixing between atomic or- 
bitals of Fe{OH)4 molecule and the ones of graphene in 
the hybrid system. In order to quantify the hybridization 
level of the occupied states, we define the quantity Cg^M : 

C./M = E<M (2) 
AO 

where C^m represents the projection of an electronic 
states into atomic orbitals (AO) belonging to graphene 
atoms (g) or to the molecule atoms (M). These states in 
the range [-3.2 ;-2.0] eV, with respect to the Fermi level, 
correspond to occupied electronic states where Cg/cu G 



[10%; 90%]. Within that interval, both systems (i.e. gra- 
phene and molecule) have a significant contribution to 
the global electronic wavefunctions which suggest a re- 
markable hybridization. This is consistent with findings 
shown in Fig. [3^, where we observe strongly polarized 
states in a region between [-3.0; -1.5] eV and also with 
the existence of additional mixed states . In particular, 
at k=0 and at -2.4 eV there is an electronic state which 
is particularly important since it represents the overlap 
of Fe 3d;^2 (13 %) and top C atom 2p^ orbital (8 %). Mo- 
reover, the contribution of this atomic orbital from the 
closest carbon atom is not present in any other mixed 
wavefunction. This indicates, for instance, the negligible 
interaction between the C atom and the oxygen atoms. 
All other Fe states are only present in pure molecular 
electronic states and not in mixed states. Conversely, we 
do not observe strongly hybridized spin down states. The 
hybridized states within [-3.2 ;-2.0] eV show mainly over- 
lap between oxygen 2px,y states and 2p2 states of car- 
bon atoms located at distances > 5.5 from Fe atom. On 
the basis of our hndings, we expect a spin up charge ex- 
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cess around the carbon atoms closer to Fe atom and an 
excess of spin down charges within a crown of carbon 
atoms around 5.5 A from Fe atom. This local magneti- 
zation within graphene plane is shown in Fig. [3]d. The 
induced magnetization in the graphene plane is sugges- 
ted by concentric crowns of spin up/down electron in 2p 
states of carbon atoms. 



IV. CONCLUSIONS 

We report a case study of an iron-based magnetic mo- 
lecule grafted on a graphene sheet and, through ab initio 
DFT calculations, we predict an induced magnetization 
effect to the graphene substrate. We also describe a typi- 



cal interaction length of 5-6 A within the graphene plane 
with respect to structural deformations, localized charge 
transfer and induced magnetization. The latter originates 
from a significant coupling between 2pz and Fe Sd;,2 or- 
bital for the closest carbon atom. On the other hand, we 
observe that the presence of an iron-based magnetic mo- 
lecule in the close vicinity of a graphene sheet does not 
strongly perturb the electronic properties of the latter, 
thus providing a promising perspective to the design of 
non-destructive sensor devices. 

We acknowledge useful discussions with X. Blase, V. 
Bouchiat, A. Candini, M. Urdampiletta and W. Werns- 
dorfer. Fig. 1 and 3 have been performed with XCrys- 
den package^'^. Calculations were done at IDRIS (Orsay, 
France), Project No. CP9-91387. 
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